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I. INTRODUCTION 



Several types of time-dependent oscillators have been studied along the past years. Ex- 
amples are: (i) the harmonic oscillator 1 ; (ii) the pseudo-harmonic oscillator 1,2 ; (iii) the para- 
metric oscillator 3 ; and (iv) the inverted harmonic oscillator 4 . Recently, another interesting 
class of time-dependent oscillators, named log-periodic oscillators, was studied 5 . 

In Ref. 5, Ozeren 5 considered the time evolution of five different one-dimensional classical 
oscillators. The coherent states for each system were constructed by using the SU(1, 1) 
algebra and their time evolution was investigated. 

In this work, we use the Lewis and Riesenfeld 6 (LR) invariant method and a unitary 
transformation to obtain the exact Schrodinger wave function for three out of the five log- 
periodic-type oscillators investigated by Ozeren 5 , namely: (i) m(t) = m -]- and k(t) = k ^f; 
(ii) m(t) = m and k(t) = k (y) 2 ; (hi) m(t) = m (j^j and k(t) = k . In all three cases 
uu(t) = Uo'f. 
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The wave functions ip n (q,t) for the time dependent harmonic oscillator (H(t) = 2 m{t) ^ 
\m{t)oj 2 {t)q 2 ) obtained in Ref. 1 are written in terms of p, a c-number quantity satisfying 
the generalized Milne-Pinney equation (p + j(t)p + oj 2 (t)p = m2 ^ p3 ), whose solution can 
be found following the procedure reported in Refs. 7 and 8. Here we write the solution 
of the Milne-Pinney equation for each system to obtain the exact wave functions for the 
oscillators. This paper is outlined as follows. In Sec. II we briefly review the LR invariant 
method for the time-dependent harmonic oscillator. In Sec. Ill we obtain the wave functions 
for the oscillators considered, and calculate the correlations between position and momentum 
and the uncertainty product. For oscillator (i) we construct the coherent states, while for 
oscillators (ii) and (iii) we construct the squeezed states. The analysis of the phase diagram 
for the three oscillators is also presented. Finally, some concluding remarks are added in 
Sec. IV. 

II. THE LEWIS AND RIESENFELD INVARIANT METHOD - WAVE 
FUNCTIONS FOR A TIME-DEPENDENT HARMONIC OSCILLATOR 

Consider a time-dependent harmonic oscillator described by the Hamiltonian 
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H(t) 



P 
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2m(t) + rW^Mo 2 . 



(1) 



whose mass {m(t)) and angular frequency (u (t)) depend on time explicitly, and the variables 
q and p are canonical coordinates with [q,p] = ih. From Eq. (1), we obtain the equation of 
motion 



where 



q + 1 {t)q + u 2 {t)q = Q ) 



l(t) = j t \nm(t). 



(2) 



(3) 



It is well known that an invariant for Eq. (1) is given by 6 



7 = 1 
2 
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where q{t) satisfies Eq. (2) and pit) satisfies the generalized Milne-Pinney 7 equation 



p + 1 (t)p + uj 2 {t) P = 



m 2 (t)p 3 



(5) 



The invariant I{t) satisfies the equation 



dl dl 1 [T TT . 



(6) 



and can be considered hermitian if we choose only the real solutions of Eq. (5). Its eigen- 
functions, (f) n (q,t), are assumed to form a complete ortho normal set with time-independent 
discrete eigenvalues, X n . Thus 



I(j)n(q,t) = \ n <f) n (q,t), 

With ((f)nAn') = 8 n n>- 

Consider the Schrodinger equation (SE) 



(7) 



(8) 



where H(t) is given by Eq. (1) with p = —ih^. Lewis and Riesenfeld 6 showed that the 
solutions vp n (q,t) of the SE (see Eq. (8)) are related to the functions (f) n (q,t) by 



^„(?,0 = e*" (t ty»(?,*), 
where the phase functions 9 n (t) satisfy the equation 



(9) 



>(?,*)>■ 



The general solution of the SE (Eq. (8)) may be written as 

n 

where c n are time-independent coefficients. 
Next, consider the unitary transformation 



(10) 



(11) 



<t>' n (q,t)=U<j> n (q,t) 



where 



(12) 



U = exp < —i 



m(t)p 
2hp 

Under this transformation and defining a = q/p, Eq. (7) now reads 



(13) 



h 2 \ d 2 



2 J da 2 



\ n (p n (cr), X n = [n+-)h, 



(14) 



where I' = Ullfi and = <p' n . The factor p 1 / 2 warrants the normalization condition 

J <l>n(q,t)<i>' n (q,t)dq = J ip* n (q,t)ip n (q,t)da = 1. (15) 

The solution of Eq. (14) corresponds to that of the time- independent harmonic oscillator 
with A n = (n + \)h . Then, by using Eqs. (12), (13) and (15) we obtain 



(/>n(q,t) 



1 



7rl/2/il/2 n |2np 



" 1/2 


im(t) / 


exp 
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(16) 



here if n is the usual Hermite polynomial of order n. 
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Applying U to the right-hand side of Eq. (10) and after some algebra, we obtain 



Finally, using Eqs. (9) and (16) the exact solution of the SE for the time-dependent 
harmonic oscillator reads 



ip n (q,t) = e 



1/2 



exp 



im{t) f p 



2h \p m(t)p 2 



x H„ 



h 



1/2 



• (18) 



III. WAVE FUNCTIONS OF TIME-DEPENDENT LOG-PERIODIC 
OSCILLATORS 



In Ref. 5, Ozeren considered five different variations of m(t) and k(t), namely: (i) 
m(t) = m and k(t) = k (y) 2 ; (ii) m(t) = m (j^j and k(t) = k ; (iii) m{t) = m (j^j 
and k{t) = k (^) (Q+2) ; (iv) m(t) = m (±) and k(t) = k (f ); and (v) m(t) = m 
and k(t) = k (j^j . Here we consider only three ((i), (ii) and (iv)) out of the five oscillators 



studied by Ozeren 5 , for which u(t) = J — u^yf. 



A. m(t) = ra ^ and k(t) = fc y 



In this case Eqs. (2) and (5) read 



•■ 1 • w o*o 



(19) 



and 



1 . uffi 



t 2 1 



P+ - p + — p= m 2 t 2 p^ 



(20) 



t r ' t 2 
respectively. 

Following the procedure described in Ref. 7, we find p = c = ^L^- . From Eqs. (17) and 
(18) we have 



^ n (q,t) = e-< n+ ^" toln ^ 



m u 



nh(n\) 2 2 



2n 



1/4 



exp 



2H 



x H n 



^m u y/ 2 ^ 



(21) 
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which, except for the phase factor, is similar to the well-known wave function for the time- 
independent harmonic oscillator. 

The coherent states for the time-dependent harmonic oscillator (Eq.(l)) are constructed 
as follows 9 . Consider the time-dependent creation (a'(£)) and annihilation (a(t)) operators 
defined as 



a t(0 



a(t) 



1/2 



2hJ 



i_y /2 

2h J 



p 



— i(pp — mpq) 



+ i(.PP ~ mpq) 



(22) 



(23) 



where [a^(t), a(t)] — 1. In terms of a(t) and at(t) the invariant I (see Eq. (4)) can be written 

as 



I = h (a\t)a{t) + ^ . 
Let \n,t) be the eigenstates of I. Therefore the following relations hold 



(24) 



a(t) = y/n\n — 1, t) 



(25) 



a\i) = Vn+ l\n+ l,t), 



(26) 



I\n, t) — H ( n + - ) \n, t). 



(27) 



Since the coherent states for I can be easily constructed, the coherent states for the 
time-dependent harmonic oscillator are straightforwardly obtained: 



oo T 
2 In & 



n=0 



(n!)V2 



i8„(t) 



n,t), 



(28) 



where 9 n is given by Eq. (17), and the complex number a(t) satisfies the eigenvalue equation 



a(t)\a, t) = a(t)\a, t), 



(29) 



with 



a(t) = a(t )e 2ie °V 



(30) 



6 



and 



W = ~ I Z^TTTy (31) 



1 /•* <ft' 

2 A, 

The fluctuations in q (Aq) and p ( Ap) and the uncertainty product ( AqAp) in the coherent 
state \a,t) , read 



(32) 



(33) 



and 



Ag a A Pa = ^(l + m 2 pV) 1/2 , (34) 

respectively. 

If m(t)pp 7^ 0, Aq a Ap a is not minimum, indicating that the coherent states \a,t) are 
not minimum-uncertainty (coherent) states. In fact, the states \a,t) for the time-dependent 
harmonic oscillator are equivalent to the well-known squeezed states, as pointed out in Refs. 
10 and 11. 

For p = c, p = and Aq a Ap a = | , indicating that the states \a,t) are "true" coherent 
states. This is an interesting result since the minimum uncertainty product is assumed to 
be satisfied only for time-independent harmonic oscillator, unless the solution of Eq.(5) is a 
constant 1 . 

Next, let us analyze the time behavior of (q) a , (p) a and the phase diagram (q) a x {p) a . 
By setting a(t ) = u + iv and using Eqs. (22) and (23), we find 



(q) a = V2h [u cos (9 Q (t)) - v sin (9 (t))] (35) 



( P ) a = V2h 



— h wrap cos (9o(t)) + ( vmp sin (9o(t)) 

P J \P 



(36) 



The constants u and v are determined from the initial conditions (q(to)) a = Qo and 
ip(to))a =Po = m(to)v . For q = 1 and v = 0, we find 



(q) a = cos \ tav In , (37) 

(p) a = -moUJo sin (t u In ^ . (38) 

Figures 1(a) and (b) show the time dependent behavior of (q) a and (p) a , respectively. 
In all plots we used t = 1.0, u = 10.0 and m = 1.0. From Fig. 1(a) we observe that 
the system oscillates forth and back between the classical turning points with an increasing 
period and constant amplitude. The phase diagram is shown in Fig. 1(c). Even though that 
this system is dissipative (total energy E — it behaves like the usual time-independent 
harmonic oscillator (E = constant). This can be seen from the relation A = y^f, where A 
is the amplitude of motion. Since k oc \ and E oc |, A is a constant. As t increases, the 
frequency u(t) decreases (oc |) while the period increases (oc ^) leading to the "exact" log 
periodic behavior shown in Fig. 1(a). 

Pedrosa et al 12 have combined linear invariants and the LR method to obtain the exact 
wave function for a particle trapped by oscillating fields, which were written in terms of 
Mathieu functions. They calculated AqAp and the quantum correlation between q and p, 
defined by C± t i = \{{qp + pq)) — (q)(p) 13 - They are related through the equation 



2 ^ 2 



AqAp=- X ll+^-C 1 , 1 j , (39) 

which shows that AqAp is minimum whenever C± tl = 0, as it happens for C^i calculated in 
the coherent state \a,t), i.e., (Ci ) i) a = 0. The fact that C ltl = does not mean that q and 
p are uncorrected. In order to verify the correlation between q and p one may study the 
function C„ i?n = \((q n p m + p m q n )) — (q n )(p m ). For the coherent state \a,t), we find that 
(£2,2)0 = — -y , indicating that q and p even assumed as "classical" quantities are correlated. 

The uncertainty product and correlations in the state ip n (Eq.(21)) are more easily cal- 
culated using the relation \ip n (q, t)) = e l9n ^\n : t). They are given by 

Aq^ n Ap^=(n + ^\h, (40) 



(41) 
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FIG. 1: Plots of (a) (q) a , (b) {p) a , and (c) the phase diagram (p) a vs (q) a . In the plots we 
used to = 1.0, q = 1.0, v = 0.0, u = 10.0 and m = 1.0. 



and 



(C 2)2 ) = -(n 2 + n+£)/l 2 . 
We noticed that Eq. (39) is satisfied for n = 0. In this case ip$ given by 



is the coordinate representation of the coherent state 14 . 



V 1/4 


mlulq 2 ' 


J exp 


2h 



(42) 



(43) 



B. m(t) = m and k(t) = k ) 2 

In this case Eqs. (2) and (5) are, respectively, given by 
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and 



(44) 



t 2 I 

2 

Following the procedure described in Refs. 7 and 8, we find p 
from Eqs. (17) and (18) we have 



(45) 

-7 ^tttt and 

-o (4^-1) 1/4 



|)(4^-l) 1/2 ln( i L) 



2+2 tNl/2 



771q (4cUgtg - 1 



1/4 



X 



f l/4 



xexp{^[,-(4^-l) 1/2 ] g 2 }x^ 



/mo \ V2 (4u#| 



- 1 



,V4 



The values of AqAp and C^i in the state |?/> n (g, £)), are given by, respectively 



(46) 



A, fc Ap fc = (4 ^ i)1/2 (n + i) ft , 



and 



(47) 



(Ci,i) t 



4u;gtg - 1 



1/2 



For n = 0, Aq A) Ap A) 



(4a,gtg-l) 



572 ft, and the state 



(48) 



^o(g,*) = e" 



i(4^-l) 1/2 l„(X 



2+2 i\ 1/2 



m (4u;gtg - 1 
2irh 



1/4 



X 



*V4 



xex P{^[ < -(^o-l) 1/2 ? 2 ]} (49) 



corresponds to the coordinate representation of the squeezed state 14 . 

For the sake of comparison with case A, let us discuss the behavior of the classical 
variables q and p on time, as well as the phase diagram. By solving Eq.(44), the solutions 
for q and p satisfying the initial conditions q — 1 and v — are, respectively, given by 



to 



, (4ujg - 1) 1/2 t 
cos in 



to (4u;gfg -if 2 



,' (4^-l) 1/2 , t 
sm, ln- 



(50) 
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FIG. 2: Plots of (a) q, (b) p, and (c) the phase diagram p vs q. In the plots we used 
to = 1-0, go = 1-0, fo = 0.0, uq = 10.0 and mo = 1.0. 



Figures 2(a) and 2(b) show the variation of q and p on time, respectively. Unlike from 
case A where the system oscillates back and forth between the turning points with constant 
amplitude, here q increases while that p decreases as time increases. Figure 2(c) shows 
the phase diagram for this oscillator. Initially at rest, the particle is speeded up, and then 
slowed down, indicating that the system is also dissipative. Since E oc 1/t and k oc l/t 2 , 
the amplitude A increases as A oc \ft. Due to the presence of the factor \ft in Eq. (50) this 
oscillator exhibits a pseudo-log-periodic behavior. 



And 




(51) 
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C. m(t) = m (j^j and k(t) = k 



Now Eqs. (2) and (5) are, respectively, given by 



2 r 
q + -q + UJ 2 ^q = 



and 



(52) 



f 4 1 



Again, by following the procedure described in Refs. 7 and 8, we find p 
and from Eqs. (17) and (18) we obtain 



(53) 
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7T^ (^!) 2 2 2n+1 



1/4 



rap 

2^o 



xt 1 / 4 

1/2 
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(54) 



The values of AqAp and Ci ; i in the state \tp n (q,t)}, are respectively given by 



^ Liu 

-if/ 2 r + 2 1A 



(55) 



and 



(C: 



i.i 



4^ " 1 



1/2 



n + - U. 



(56) 



The expression of the coordinate representation of the squeezed state for n = 0, reads 



,1/2 
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■ 111 
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(57) 



By solving Eq. (52) and using the initial conditions q — 1 and t> = 0, we find 
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FIG. 3: Plots of (a) q, (b) p, and (c) the phase diagram p vs q. In the plots we used 
to = 1-0, go = 1-0, t>o = 0.0, uo = 10.0 and m = 1.0. 

and 

V*o(4^-l) 1/2 V 2 *o 7 

The behavior of the classical g and p variables on time is displayed in Figs. 3 (a) and 
(b), respectively. Despite the oscillating (cos ^^lntj and (sin ^^lntj) terms, q and p 
exhibit an opposite behavior compared to those calculated in case B. Here q decreases 
while p increases as t increases. Figure 3(c) shows the phase diagram. We observe that 
the amplitude A decreases as oc l/y/i. This oscillator also exhibits a pseudo-log-periodic 
character. 
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IV. CONCLUDING REMARKS 



In this paper we have used a unitary transformation and the LR invariant method in 
the Schrodinger picture to obtain the exact wave functions for oscillators exhibiting either 
log-periodic or pseudo-log-periodic behavior. It is well-known that a challenge in obtaining 
the exact solution (see Eq. (18)) for the SE (Eq.(8)) for H(t) given in Eq. (1), is the solution 
of the auxiliary equation for the c- number quantity p (see Eq. (5)). Here we find p for each 
case using the methods described in Refs. 7 and 8. 

For case A, we find p = c and, as a consequence, the solution for ip n (q,t) (see Eq. (21)) 
except for the phase factor ( e ~^ n+ 2)* (,Woln ^) j is similar to the well-known wave function 
for the time-independent harmonic oscillator of mass m and frequency uq. In Ref. 1, we 
observed that when m(t) = m , u(t) = u , and p(t) = ( -^j-J , which is a particular solu- 
tion of Eq. (5), the wave function obtained also corresponds to that of the time-independent 
harmonic oscillator. In case A even with m oc t and to oc |, we obtain the same solution for 
p (p = c), indicating that this oscillator behaves as the harmonic oscillator with m and u 
constant. 

We have constructed the "true" coherent states, \cx,t), whose coordinate representation 
is given by Eq. (43). We verified that Eq. (39) holds for \a,t). We calculated the quantum 
fluctuations in the coordinate and momentum as well as the quantum correlations between 
the coordinate and momentum in the state ijj n {q,t). 

We analyzed the time behavior of (q) a and {p) a , as well as the phase diagram (q) a x (p) a 
(see Fig.l(a)-(c)). We observed that (q) a and (p) a exhibits the exact log-periodic behavior, 
and that the phase diagram indicates, as already anticipated, that the log-periodic oscillator 
behaves as the classical harmonic oscillator with m(t) = m and ou(t) = u . 

For cases B and C, we obtained the wave functions given by Eqs. (46) and (54), respec- 
tively. 
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